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The  variance  <r|  of  the  Strehl  ratio  of  a  reasonably  well-corrected  adaptive  optics  system  is  derived  as  a  power 
series  in  the  log-amplitude  variance  (rf  and  the  residual  phase  error  variance  ( r It  is  shown  that,  to  leading 
order,  the  variance  of  the  Strehl  ratio  is  dependent  on  the  first  power  of  the  log-amplitude  variance,  (of)1,  of 
the  incident  optical  field  but  only  on  the  second  power  of  the  residual  phase  variance,  (cr^) 2 ,  of  that  field  after 
adaptive  optics  correction,  and  on  the  first  power  of  the  product  of  the  log-amplitude  variance  times  the  phase 
variance,  ( (rf<r%h) 1 .  As  long  as  the  adaptive  optics  correction  is  good  enough  to  ensure  that  the  variance  of  the 
residual  phase,  is  significantly  less  than  unity,  then  even  for  fairly  small  values  of  the  log-amplitude  vari¬ 
ance  of,  the  value  of  the  variance  of  the  Strehl  ratio,  of,  will  be  dominated  by  the  value  of  the  log-amplitude 
variance.  ©  1998  Optical  Society  of  America  [S0740-3232(98)02008-0l 
OCIS  codes:  010.1080,  010.1330,  290.5930,  350.5030. 


In  studying  the  variability  of  the  Strehl  ratio  of  an  adap¬ 
tive  optics  system,  we  have  derived  an  analytic  expression 
for  the  variance  of  the  Strehl  ratio,  rr|,  as  a  power  series 
in  the  residual  phase  variance  and  the  log-amplitude 
variance  <rf  in  a  regime  where  each  of  these  two  variances 
is  significantly  less  than  unity.  We  have  been  surprised 
to  find  that  our  result  for  <r%  has  no  linear  dependence  on 
the  residual  phase  error  variance  of  (/j.  The  consequence 
of  this  is  that  in  general  the  value  of  the  variance  ct%  of 
the  Strehl  ratio  is  dominated  by  the  value  of  the  incident 
optical  field’s  log-amplitude  variance  of  and  not  by  the 
value  of v,  of  the  variance  (or  mean  square  value)  of  the 
residual  phase  error  after  adaptive  optics  correction. 
While  improving  the  adaptive  optics  performance,  i.e.,  re¬ 
ducing  the  value  of  of  f/j,  will  improve  the  mean  value  of 
the  Strehl  ratio,  (S),  bringing  it  closer  to  unity,  beyond 
some  fairly  modest  level  of  performance  the  further  im¬ 
provement  of  the  adaptive  optics  will  not  significantly  re¬ 
duce  the  variance  of  the  Strehl  ratio.  At  that  point  the 
value  of  the  variance  of  the  Strehl  ratio  is  dominated  by 
the  pattern  over  the  system’s  aperture  of  the  incident 
field’s  optical  intensity  variations— by  what  we  call  the 
“random  apodization”  of  the  aperture. 

We  have  considered  a  ground-based  adaptive  optics 
telescope  system,  a  system  looking  upward  and  operating 
in  the  visible  or  the  near  infrared— for  which  case  the  log- 
amplitude  variance  will  generally  be  of  the  order  of  of 
-  0.01  to  0.10  Np2  —  and  have  concerned  ourselves  with 
the  variability  from  instant  to  instant  of  the  system’s  per¬ 
formance,  the  performance  being  measured  by  the  Strehl 
ratio  S.  We  have  restricted  our  attention  to  the  Strehl 
ratio,  as  this  quantity  serves  as  a  surrogate  for  the  less 
exactly  defined  concept  of  resolution1  and  as  a  surrogate 


for  the  concept  of  antenna  gain.  Our  interest  in  the  ran¬ 
dom  variability  from  instant  to  instant  has  to  do  with  ap¬ 
plications,  such  as  short-exposure  imagery,  for  which  its 
value  at  a  particular  instant  that  will  determine  the  sys¬ 
tem  user’s  assessment  of  the  system’s  performance. 

The  value  of  S  is  set  by  both  the  residual  phase  error, 
i.e.,  by  the  (small)  part  of  the  turbulence-induced  wave- 
front  distortion  that  for  various  practical  reasons  the 
adaptive  optics  did  not  correct,  and  by  the  variations  over 
the  telescope’s  aperture  of  the  intensity  of  the  received  op¬ 
tical  signal.2 

An  analysis  of  the  expected  values  of  the  first  and  sec¬ 
ond  moments  of  the  Strehl  ratio,  (S)  and  (S2),  respec¬ 
tively,  shows  that  these  quantities  have  values  given  by 
the  equations 

<S> 


JJ  dr,  ,dr2 


*2W(r]  )W(r2)exp(Fj) 


(i) 


dr2W(r,)W(r2 


<S2> 


////  dr]  ,dr2dr3dr4W'(r1  )W(r2) W(r3)W(r4)exp(F2) 
j  j  j  j  d^  ,dr2dr3dr4W(r1)Wr(r2)W(r3)W(r4) 

(2) 

Although  these  equations  were  derived  before  the  devel¬ 
opment  of  adaptive  optics3  and  the  wave-front  distortion 
considered  then  was  the  uncorrected  wave-front  distor¬ 
tion  (i.e.,  wave-front  distortion  without  any  adaptive  op- 


2108 


H.  T.  Yura  and  D.  L.  Fried 


tics  correction),  those  results  are  fully  applicable  to  an 
adaptive  optics  system  if  the  measure  of  the  wave-front- 
distortion  statistics  inherent  in  the  Fj  and  F2  functions 
appearing  in  Eqs.  (1)  and  (2)  is  understood  to  be  that  for 
the  residual  wave-front  distortion  after  adaptive  optics 
correction. 

The  function  W(r)  denotes  the  extent  of  the  system’s 
aperture,  having  a  value  of  unity  for  a  position  r  that  is 
inside  the  aperture  and  a  value  of  zero  for  a  position  that 
is  outside  the  aperture.  We  have  considered  the  system 
that  has  a  clear  circular  aperture  of  diameter  D ,  for 
which  case  we  can  write 

fl  if  |r|  *£  D/2 

W<r)  =  (o  if  |r|  >  D/2  '  (3) 

The  arguments  of  the  exponential  functions  in  Eqs.  (1) 
and  (2),  denoted  by  Fx  and  F2,  have  values  given  by  the 
equations 

Fi  =  -kl'Ssd rj  -  r2)  +  './/(r,  -  r2)]  (4) 

F2  =  -k['Ss*iT\  -  r2)  -  '•/**(*  1  -  r3)  +  '-/*<*!  -  r4) 

+  '/**(  r2  -  r3)  -  </S(h{  r2  -  r4) 

+  '/>*(r3  -  r4)  +  '//(r!  -  r2)  -  '//(  rj  -  r3) 

+  '_//( r,  -  r4)  +  y;(r2  -  r3)  -  f/,(r2  -  r4) 

+  .//( r3  -  r4)]  +  2[‘/,(r1  -  r3)  +  ^/(r2  -  r4)]. 

(5) 

The  notations  r)  and  V/(r)  denote  the  structure 
functions  for  the  (residual)  phase  and  for  the  log  ampli¬ 
tude.  The  notations  ‘^^(r)  and  5</(r)  denote  the  covari¬ 
ance  of  the  residual  phase  and  of  the  log  amplitude,  re¬ 
spectively. 

The  structure  functions  and  the  covariance  functions 
are  related  by  the  equations  r)  =  2[^^(0) 

"  ‘^^(r)]  and  l//(r)  =  2[‘//(0)  -  K/( r)],  and  the  vari¬ 
ances  and  the  covariance  functions  are  related  by  the 
equations  cr2Sth  —  ‘^^(0)  and  of  =  ^/(  0).  We  can  ex¬ 
press  the  structure  functions  as 

-/<5</»(r)  =  2o^[l  —/***<  r)], 

where  /^(r)  =  (6) 

and 

'//(r)  =  2of[l  -  //( r)],  where  /)( r)  -  ^/(r)/of. 

(7) 

With  this  notation  it  is  convenient  to  recast  Eqs.  (4)  and 
(5)  as 

F{  =  +  fs+i* i  ~  r2)]  +  "f[-l  +  //( rx  -  r2)], 

(8) 

F2  =  <4j‘2  +  /^(ri  -  r2)  -  /^(r,  -  r3) 

+  -  r4)  +  r2  -  r3)  -  /^(r2  -  r4) 

+  /^<*3  “  r4)]  +  "f[“2  +  fi( r,  -  r2) 

+  //( r j  -  r3)  +  f,(r]  -  r4)  +  />( r2  -  r3) 

+  /)(r2  -  r4)  +  //( r:J  -  r4)l.  (9) 


Based  on  the  assumption  that  the  variance  of  the  residual 
phase,  of  f/„  and  the  log-amplitude  variance  of  are  both 
small,  and  recognizing  that  /^(r)  and  fj( r)  have  magni¬ 
tudes  no  greater  than  unity,  we  consider  the  power  series 
expansion  of  the  exponential  functions  in  Eqs.  (1)  and  (2), 
retaining  only  terms  up  to  second  order  in  the  argument 
of  the  exponential,  i.e.,  using  the  approximation  that 
exp(x)  1  +  x  +  \x 2.  With  this  approximation  we  ob¬ 
tain  the  results  that 

exp(F, )  ~  1  +  { (t%,[  ~  1  +  r,  -  r2)] 

+  of[-l  +  fi(r1  -  r2 )]} 

+  -  %fs4ri  ~  r2> 

+  [fs,i,(r i  -  r2)]2}  +  (of)2{  1  -  2/‘,(r1  -  r2) 
+  [ //( ri  ~  r2)]2}  +  2(t%<t21[1  -  fs<ll(ri  -  r2) 
-  //( rj  -  r2)  +  fH( rj  -  r2)f,(r1  -  r2)]) 

(10) 

and  (after  considerable  algebraic  manipulation)  the  result 
that 

exp(F2)  =>1  +  {a%[-2  +  2f^(rx  -  r2)] 

+  erf [-2  +  6f,( r!  -  r2)]} 

+  l{(o-w>2[4  -  8fs4 rj  -  r2)  +  Kw] 
+  (of)2[4  -  24f/(r1  -  r2)  +  K{] 

+  2<4/,of[4  -  4fSll>(  r ,  -  r2) 

-  12 /-/(r,  -  r2)  +  (11) 

where 

Ku  =  6/Vr,  -  r2)fsJ,(r1  -  r2) 

-  8fsJ *•]  -  r2)fs,h(r2  -  r3) 

+  6/'w( r,  -  r2 )/'rfu!,( r3  -  r4),  (12a) 

K,  =  6//( r,  -  r2)fi(r1  -  r2)  +  24/',(r1  -  r2)f,(r2  -  r3) 
+  6/)(r,  -  r2)fi( r3  -  r4),  (12b) 

K,UJ  =  2  fu(  r,  -  r2  )/>(**,  -  r2)  +  8/'^(rJ  -  r2) 

X  //(r2  -  r3)  +  2 /",)•,/,( r,  -  r2)f,(r3  -  r4).  (12c) 

When  relations  (10)  and  (11)  are  substituted  into  Eqs.  (1) 
and  (2),  we  can  obtain  the  results  that 

(S)  =  1  +  <r2,A(-l  +  au)  +  arf(- 1  +  a,) 

+  5('rl/>)2<1  _  2«,w,  +  f},u,) 

+  |(<r2)2(  1  -  2a/  +  /3,) 

+  "■2,/,"72<  1  “  ~  “/  +  Ps.th  (13) 


and 
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(S2)  =  1  +  2  +  2a**)  +  ( —2  +  6a/) 

+  ( "‘5<*)2< 2  -  4aw  +  30m  ~  4r^  + 

+  (<t,2)2(2  -  12a,  +  3 /?,  +  12y/  +  3«f) 
+  <w,°f<4  ~  4aw  -  12a,  +  2/3,* , 


+  8  7*,/  +  2a**a,), 


where 


1dr2W(rJ)W(r2)/’,(rj  -  r2) 


J  J  dr,  dr2  W(  r, )  W(  r2 )/^( r ,  -  r2) 

JJ  dr,dr2)V(rj)W(r2) 

J  J  dr1dr2W(rJ)W(r2)/’,(r1  -  r2) 

J  J  dr,dr2W(r,)W(r2) 

J  J  dr]dr2W(ri)W(r2)[/'^(rI  -  r2)]2 

JJ  dr,dr2W(r,)W(r2) 


SS dri 


,dr2W(r,)W(r2) 


dr2W(r,)W(r2)[/-,(r,  -  r2)]2 


JJ  dr, 


dr2W(r,)W(r2) 


dr,  dr2W(  r,  )W(r2  )/■,,*(  r,  -  r2)/>( r,  -  r2) 


8d>,l  ~ 


//dr, 


dr2W(r1)W(r2) 


The  variance  of  the  Strehl  ratio,  <r|,  is  expressed  in  terms 
of  the  first  and  second  moments  of  the  Strehl  ratio,  (S) 
and  (S2),  respectively,  by  the  equation  cr|  =  ( S 2) 
-  <S)2.  Making  use  of  Eq.  (13)  and  after  carrying  out 
the  squaring  process  and  dropping  terms  of  too-high  or¬ 
der,  we  obtain  the  result  that 

(S)2  =  1  +  (T%(  -2  +  2 aH)  +  of(- 2  +  2a,) 

+  (a^)2(2  -  4 au  +  a%  +  pH) 

+  ( erf )2(2  -  4a,  +  af  +  /?,) 


+  < T2s4,<rf(4  -  4 aS4, 


4a/  +  2a^a/  +  2/?**/). 


This  result  together  with  Eq.  (14)  allows  us  to  write  for 
ogj  the  variance  of  the  Strehl  ratio,  that 

rr|  =  of(4a/)  +  (o**)2(2/3**  +  2 a2S(f)  -  4 y<*) 

+  (<x2)2(  —  8a/  +  2/6/  +  2a2  +  12y/) 

+  ^<7/2(_8a/  +  (19) 

Considering  this  result,  we  can  see  that  the  variance  of 
the  Strehl  ratio  is  first  order  in  the  log-amplitude  vari¬ 
ance  but  only  second  order  in  the  variance  of  the  residual 
phase.  We  can  consider  this  our  final  result,  or  if  we 
wish,  in  the  spirit  of  order-of-magnitude  approximations, 
we  can  make  the  approximations  that 

-  ( d/D )2,  a,  -  (p/D)7'3,  (20a) 

P„~ld/D)*,  ~  (p/D)7/3,  0MJ~  (d/D)2,  (20b) 


yM  ~  (d/D)4,  7/  -  (p/D)14/3,  7*,/  ^  ( d/Df(P/D )7/3, 

(20c) 

where  cf  denotes  the  adaptive  optics  subaperture  size  and 
may  be  considered  to  be  the  correlation  length  for  the  re¬ 
sidual  phase,  and  p  denotes  the  correlation  length  for  the 
log-amplitude  variations.  These  approximations  follow 
from  consideration  of  Eqs.  (6)  and  (7)  and  Eqs.  ( 15)— (17) 


Jsihj  ~ 


J  J  J  dr,  dr2dr3W(  r1)W7(r2)W(r3  )/"**(  rj  -  r2  )fM(  r2  -  r3) 
JJJ  dr,dr2dr3W(r,  )W(r2)W(r3) 

JJJ  dr, dr2dr3 W( r , )  W(  r2 )  W(  r3 )/,( r,  -  r2)f,(r2  -  r3) 
JJJ  dr,dr2dr3W(r,)W(r2)W(r3) 

JJJ  drldr2dr3W(rI)W(r2)W(r3)fM(r1  -  r2)f,( r2  -  r3) 

JJJ  dr,dr2dr3W(r,)W(r2)W(r3) 


(17a) 


,  (17b) 
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and  from  consideration  of  the  facts  that  the  value  offS(h(r) 
goes  to  zero  rather  rapidly  for  values  of  r  for  which  |r| 
>  </,  that  the  value  of  /}(  r)  goes  to  zero  (some  what  more 
slowly)  when  the  value  of  r  is  such  that  |r|  >  p,  that  for 
all  adaptive  optics  systems  d  D ,  and  that  the  value  of  p 
is  of  the  order  of  a  few  to  a  few  tens  of  centimeters  (de¬ 
pending  on  the  wavelength),  while  the  telescope  aperture 
diameter  D  is  of  the  order  of  meters— so  p  <§  D.  That 
the  p  over  D  ratio  shows  up  in  relations  (20)  raised  to  the 
seven-thirds  power  rather  than  to  the  second  (or  six- 
thirds)  power  is  a  consequence  of  the  way  in  which  /}( r) 
approaches  zero  as  the  value  of  |r|  increases— with  an  ini¬ 
tial  undershoot  (i.e.,  a  drop  below  zero  value)— as  devel¬ 
oped  by  Gracheva  and  Gurvich4  and  by  Lutomirski  et  a/.J 
In  writing  the  expression  for  ps<hj  we  have  made  use  of 
the  fact  that  (except  under  conditions  of  saturated 
scintillation6)  p  will  be  greater  than  r0  and  thus  greater 
than  d,  so  the  effective  range  of  the  integrand  in  Eq.  (16c) 
is  set  by  d  and  not  by  p. 

The  approximation  sign  is  used  in  relations  (20)  to  in¬ 
dicate  only  an  order-of-magnitude  dependence.  In  this 
order-of-magnitude  sense  we  can  recast  Eq.  (19)  as 

(rs  **  p/Z) ) d/D  )'2  +  ( (rf  )2{p/D)'1/,i 

+  (21) 


just  as  we  might  have  expected.  The  leading  value  in  the 
deviation  from  unity  of  the  mean  value  of  the  Strehl  ratio 
can  be  either  (r  %t  or  a  f  depending  on  which  is  greater. 
The  deviation  from  unity  of  the  mean  value  of  the  Strehl 
ratio  can  be  dominated  by  the  mean  square  residual 
phase  error,  whereas  the  variance  of  the  Strehl  ratio  can¬ 
not. 

From  consideration  of  relation  (23)  we  can  see  that  ( S ), 
the  mean  value  of  the  Strehl  ratio,  is  dependent  on  the 
first  power  of  a\h,  the  mean  square  value  of  the  residual 
phase  error.  The  effect  of  of,  the  log-amplitude  vari¬ 
ance,  which  generally  is  significantly  smaller  than  a\h,  is 
consequently  generally  only  of  secondary  importance  in 
determining  the  value  of  (S).  On  the  other  hand,  we  see 
from  consideration  of  Eqs.  (21)  and  (22)  that  for  o-f,  the 
variance  of  the  Strehl  ratio,  generally  the  dependence  is 
principally  on  of,  the  log-amplitude  variance,  with  only  a 
second-order  dependence  on  of5(/;,  the  mean  square  value 
of  the  residual  phase  error— even  though  the  value  of  a%} 
can  be  significantly  larger  than  the  value  of  of. 
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The  first  thing  to  be  noted  in  considering  Eq.  (21)  is  that 
in  all  existing  adaptive  optics  systems  d  <  D  and  that  for 
all  propagation  conditions  of  interest  the  value  of  p  is 
comparable  to  or  larger  than  r()  and  is  thus  comparable  to 
or  larger  than  d . 

This  allows  us  to  conclude  that  to  leading  order,  the 
Strehl  ratio  variance,  erf  does  not  depend  on  the  variance 
of  the  residual  phase,  (r%f),  but  rather  on  the  log- 
amplitude  variance  of.  Only  in  the  extreme  case  where 
the  scintillation  is  so  weak  that  of  is  less  than 
( of(/j)2( d/p)2  would  the  (  of (/> ^-dependent  term  be  domi¬ 
nant.  The  general  form  of  the  order-of-magnitude  result 
can  be  written  as 


<rf{p/D)1:i  if  of  ^  itrih)z{dlp)HDfp)m 

( aiAd/D  )2  otherwise 

(22) 


This  presumes  for  its  development  that  ofv,  <  1. 

Starting  from  Eq.  (13),  the  corresponding  results  for 
the  mean  value  of  the  Strehl  ratio,  (S),  is  given  by  the  re¬ 
lation 

<s>  -  1  -  <4,  -  '4  (23) 
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